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Abstract 



The proton structure effects, including finite size, polarizability and self- 
Q_i| energy is considered and their influence on energy levels of muonic hydrogen 

is recalculated. A new theoretical prediction for the Lamb shift is presented 
together with improved values of all known QED contributions. 
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The precision tests of quantum electrodynamics (QED) in atomic systems have reached 
the level, where nuclear structure effects become significant. Moreover, the lack of accurate 
data on low energy structure functions of the nucleus strongly limits theoretical predictions. 
The well known example is the hydrogen hyperfine structure splitting, where measurement is 
6 orders of magnitude more precise than the current theoretical predictions. Other example 
is the Lamb shift in hydrogen, where inaccuracy in the proton charge radius dominate other 
theoretical uncertainties. Since many years it was the motivation to study the pure QED 
systems like muonium and positronium, where strong interaction effects are negligible or well 
estimated at the precision level of interest. A new possibilities in the improvement of QED 
test on bound systems appeared with the project of measurement of the Lamb shift in muonic 
hydrogen Jl|. Being bound in the ground state, muon with its 200 heavier mass compared 
to the electron, penetrates the proton and become sensitive to the distribution of charge, 
magnetic moment or the polarizability. The main goal of the measurement of the Lamb 
shift in muonic hydrogen is a precise determination of mean square proton charge radius. 
It will verify older measurements of proton charge radius based on the elastic electron- 
proton scattering. The results of these experiments had to be extrapolated to very low 
q 2 = momentum transfer, and therefore might be not very reliable. In contrast the 
muonic Lamb shift measurement aims to improve the precision of the proton charge radius 
by at least 10 times. However it would require the calculation of all other contributions 
with the comparable or better precision. Few years ago, in the summary of all known up 
to date QED effects on muonic hydrogen energy levels [fj], we concluded that three-loop 
vacuum polarization have limited the precision of theoretical predictions by 0.01 meV, and 
pointed attention to the proton polarizability effect, which was not calculated that time. 
The progress on the experimental side with the preparation of the measurement of muonic 
hydrogen Lamb shift has stimulated further theoretical works. The rather difficult three- 
loop vacuum polarization contribution has recently been calculated by Kinoshita and Nio 
H with the result 

AE(2P-2S) =0.0076 meV. (1) 

A new contribution coming from hadronic vacuum polarization has been studied by Friar 
and collaborators Q, who obtained the following result 

AE(2P -2S) = 0.011 3(3) meV . (2) 

The proton polarizability correction has been analyzed by Rosenfelder in ||. His estimate 
for this effect is 

AE(2P - 25) = 0.017(4) meV . (3) 

Another recent work || gives a similar result of ~ 0.018 meV. The calculation of proton 
polarizability effect is affected by the lack of precise data on low-energy proton structure 
functions. In this paper we present another estimate of this effect together with the complete 
review of all other proton structure contributions. We would like to emphasize the impor- 
tance of the proton self-energy effect and the related problem with the meaning of mean 
square charge radius. In other words we present the relation of charge radius as obtained 
from atomic spectroscopy measurements with that based on the electron scattering data. 



It is well known, the shift of atomic energy levels due to finite charge distribution of 
atomic nucleus. It is given by the formula 

E FS = ^a 4 fi 3 (r 2 )5 l0 , (4) 

which gives —3.862(108) meV contribution to 2P-2S splitting, were we used r = 0.862(12) 
fm from ]7| . From this one concludes, that the measurement of muonic hydrogen with the 
precision of 0.01 meV will lead to the tenfold improvement in the proton charge radius. This 
formula in (^) accounts for most of the proton structure effects. Any correction beyond that 
is much smaller and usually neglected for light (electronic) atoms. It is our aim to review 
them in the context of muonic hydrogen Lamb shift. 

Since the leading contribution to the 2P-2S splitting comes from the electron vacuum 
polarization (e.v.p.), the second order corrections from combined finite size and e.v.p. effects 
are nonnegligible. They have been calculated in ||, here we present a little more accurate 
result 

AE(2P - 2S) = -r 2 0.0282 = -0.0209(6) meV . (5) 

Further corrections are due to pure photon exchange terms only. The 0(ma 5 ) correction 
is given by two-photon scattering amplitude with external momenta on mass shell, see Eq. 
(0). It is the main contribution, which we analyze in this work. However, lets consider 
first the small corrections beyond this two-photon exchange approximation, namely that 
of order ma 6 . Since they are small, it is sufficient to calculate them in the external field 
approximation, or in other words in the limit of infinite nucleus mass. The energy levels are 
obtained then from the Dirac equation in the potential of the finite size nucleus. They have 
been considered in details by Friar in |J. In the logarithmic approximation they are given 
by two terms 

£ = £ FS a 2 lna(J/i 2 (r 2 >-l) (6) 

and contribute the amount of —0.000 9 meV to the 2P-2S splitting and thus are almost 
negligible. 

The correction given by the two-photon exchange is of the main interest. It is given by 
the following expression 
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t = (1, 0, 0, 0), P = Mt is a proton momentum at rest and v = q°. For a point-like proton 
TV = t^(M) and 
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Since the amplitude T^ does not have any singularities at small q 2 the following holds: 
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The off-shell (spin averaged) forward Compton amplitude T^" of the proton is not directly 
measured. However, it could be expressed in terms of its imaginary part, through the 
dispersion relations. Before using them, one notices that the integral in Eq. (|7|) is infrared 
divergent. It requires subtraction of the leading finite size term, which has already been 
accounted for. It also requires an additional subtraction of the proton self-energy term. 
Since it changes the analytic behavior of T% and T 2 at small p 2 — M 2 , see Eq. (EDI), this 
proton self-energy could be only partially accounted for in the proton formfactors, as a 
contribution to anomalous magnetic moment or the charge radius. Since, in our opinion 
this problem is not generally known, we repeat here the analysis from our former paper on 
radiative recoil corrections and correct some minor missprints. 

If we assume a point-like proton the contribution of the proton self-energy to the Lamb 
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For a true proton there is a finite size correction as given by Eq. (|4]). The problem is that 
the proton self-energy is modified by and modifies as well, the finite size effect. Therefore 
some corrections might be counted twice. To incorporate the correction ( |i~4"| ) unambiguously 
we must precisely specify what is the nuclear mean square charge radius. Its usual definition 
through the Sachs formfactor 



(r 2 ) dG E (q 2 
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is not correct at our precision level, because the radiative correction to Ge is infrared 
divergent or depend on spurious photon mass. Following 0] we propose thus a different 
definition using the forward scattering amplitude described by T Ml/ , or more precisely by its 
longitudinal component T L defined by 



Tr 
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For our purpose we consider a nonrelativistic limit v ~ q 2 and p 2 — M 2 = (P + q) 2 — M 2 p» 
of Tl. For a point-like particle without radiative corrections Tl is 
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where G^ = F\ + -A^, F2 is an electric formfactor. The radiative corrections for a point-like 
particle || in the nonrelativistic limit are 

A ^ = 7^(y + 3 ln M^) + ( ^- 9) - I 2 ") 

We define (r 2 ) by the following equation that describes the low-energy behavior of the 
forward scattering amplitude 

m g 2 /4a , M 2 2 lrf/ ,.\ . 

We expect that for any nucleus the logarithmic term above will be the same, since it is only 
related to the fact that nucleus has a charge, and does not depend on other details like the 
spin. The associated correction to the energy for S-states has the form 
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The small second term in the above Eq. gives —0.0099 meV for 2S state in /iH. For P-states 
the proton self-energy contributes only through Bethe log. The correction to anomalous 
magnetic moment is already taken into account in the calculation of relativistic effects as 
given by the Breit Hamiltonian. 

We now return to the forward scattering amplitude and its associated correction to energy 
of //H. We can neglects here QED effects on the structure functions, since they are a times 
smaller. Otherwise dispersion relations for T\ and T 2 as in Eq. (p9jj30"D will not be correct. 
After neglection of QED effects one can expect a separated pole at 2Mu = q 2 , that is due 
to the elastic contribution. It is obtained with the help of the proton elastic formfactors F\ 
and F 2 , see Eq. flUf ). One derives for T\ and T 2 the following expressions 

T ° = ~ („4 _ l M - 2 -T) < 4 MVF l + 2 1 1 *i F 2 + q l F|) (23) 

t?= «-Im>*) ^ **-**$■ (24) 



Any contribution in the two photon exchange beyond this second order Born term is con- 
sidered to be due to the proton polarizability. The correction to energy as given by Eq. 
(0) shows infrared singularity, which is due to the finite size effect, already accounted for at 
the lower order in a. One first subtracts this singularity and then integrate over q. With 
the commonly used dipole parametrization of proton formfactors and A 2 = 0.71 GeV 2 one 
obtains AE = 0.018 meV. However, as pointed out in |l[ this result is not correct, because 
this dipole parametrization corresponds to the proton radius r = 0.81 fm, which differs 
significantly from the value r = 0.862 fm. Instead, we use the parametrization by Simon et 
at. with the result 

AE = 0.0232(15) meV = r 3 0.0281. (25) 

This parametrization in [7j is not correct for Q 2 = —q 2 « 1 GeV or greater since for large 
Q 2 it behaves like Q~ 2 not Q~ A , however it does not influence the result at the presented 
precision level. This r 3 dependence in Eq. (]25 ) is only an approximate dependence. This 



dependence become exact in the large nucleus mass limit with the dipole parametrization 
of proton formfactors. The result 0.0232(15) was obtained assuming r = 0.862(12) fm. If 
the proton radius r from /xH measurement will be significantly different, this result should 
be adjusted with the improved parametrization of proton formfactors. 

Let us assume now, that the second order Born contribution is subtracted from T\ and 
T 2 . The proton polarizability correction Epol as given by Eq. (|7|) is split into two parts 

Epol = E ± + E 2 , (26) 

E l = -2*mZ[ *L^ + ?-,.M S), (27) 



M J (2 7r) 4 i <f(g 4 -4m 2 z/ 

. 9. . m f d A q q 2 + 2 v 2 / . 9 . v 2 m . 9 , . 

E 2 = 2e 4 2 (O) T7 / 77-^- -tA 7-^7^ (T 1 (v,q 2 ) + -T 2 (u,q 2 ) . (28) 



M J (2vr) 4 i g 4 (g 4 -4m 2 z/ 2 ) V v ' q 

The first dominant part could be well calculated, while the calculation of the smaller second 
part will require additional assumptions. This second part was usually neglected in studies 
of polarizability effect in composite nuclei, see for example ||11|| . The imaginary part of Tj 
in variable v at fixed q 2 is obtained from the inclusive cross section 7* + p — > X. The real 
part could be restored from the imaginary one using the following dispersion relations: 



T 2 {v,q) = - „ du — — , (29) 
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where v t h is the threshold value of v for the production of tc mesons. This formulas are 
obtained from [p^] , by multiplying right hand side by —1/2, due to the different definition 
of Ti in our work. While W 2 for proton behaves like 1/v to large v, W\ goes like u, and thus 
the dispersion relation in (|30|) contains a subtraction at v = 0. We have assumed here, that 
T 2 has similar asymptotic behavior as W 2 . i.e. vanishes in the large v limit. However, there 
is no too much information on Ti(0,g 2 ). One recognizes that in the limit of small q 2 it is 
given by the magnetic polarizability (3m 
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which amounts to (3m = 1.56(57) • 10~ 4 fm 3 , see |13|]. Since there is no experimental data, 
we assume here, that the g- dependence is governed by the square of the elastic formfactor, 
namely 



Pm(Q 2 ) = Pm 
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(A 2 + Q 2 ) 
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with A 2 = 0.71 GeV 2 . Using equations (|27| - |3~TD one derives the following expressions for E\ 
and E 2 . 
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and R is equal to R — Wl{v, q 2 )/Wi(i', q 2 ) = (Tl/<7t, the ratio of longitudinal to transverse 
cross sections, and thus could be measured. A recent work [14]| by E143 collaborations brings 



the most recent parametrization of R in the momentum range Q > 1 GeV. What we find, 
is a weak dependence of R on v for small Q 2 , of order few GeV, see for example Figure 3 of 
this work. Moreover R is between 0.3 — 0.4. However there is no precise data on the most 
important region for us of Q 2 < 1 GeV. One knows from conditions (|T2| , [T3"| ) , that R 
small Q 2 and constant v. The only known paper by Brasse et. al 



Q 2 at 

15| brings comments on 



R, that it is of order 10% at the resonance region. Since, what we find later on, E 2 is small 
and we do not have precise data on R at the resonance region we simply neglect R in the 
expression for E 2 in (0). We expect, this approximation will not alter significantly the final 
result for the polarizability correction. Various parametrization of structure functions W 2 
were presented in the literature. We have chosen that, which match the q 2 = limit known 
as a photoproduction and are correct for region of v and Q 2 of order few GeV. We use the 
parametrization from [T5| for the resonance region. It is given by a compact formula and a 
long table of numbers, see Table 1 of |15]. For v close to the threshold for 7* + p + — > 7r + + n 
process, we assume a behavior \fW — W t h and match the value at W = 1.110 GeV, as given 



in |15| . The photoproduction of n° could be neglected close to this threshold value. Out of 



the resonance region, W > 1.990 GeV, we use an ALLM97 parametrization with the recent 
corrections |I6 |. With these parametrizations integrals in ( |33"| , |3~4T) are performed numerically 



with results 

Epol = E l + E 2 = -a 2 ^Jt (1.82 + 0.12 - 0.26) = -a 2 ^-^1.68(33) , (38) 

mm mm 

where the uncertainty of 0.33 forms 20% of final result and is an estimate due to approxi- 
mations and assumptions performed during this derivation, namely: unknown R, unknown 
Q 2 dependence of (3m and inaccuracy of ALLM97 parametrization at low Q 2 and v. The 
contribution to the (2P-2S) splitting of n = 2 states is 

E POL = 0.012(2) meV. (39) 

It is slightly lower than the result obtained by Rosenfelder in ||. We think it is due to 
approximate treatment of the Q-dependence in his work. Other contributions to the Lamb 
shift in muonic hydrogen are presented in the Table I. It is an improved version of the former 
Table I in || . All QED corrections are recalculated with a better precision. The final result 
for the Lamb shift differs from that in M by more than the error estimate there. It is due 
to the mistake in the sum in Table I in || and due to the inclusion of new terms: proton 
polarizability, third order electron v. p. and the hadronic v. p. It should be noted here that 
the same effects should be included in hydrogenic Lamb shift, when having a new value 
of the proton charge radius. It was not necessary so far since in the electron scattering 
measurements of proton radius, these effects have not been excluded. The uncertainty in 
the final result has two main sources: the proton polarizability and the estimate of higher 
order QED effects. While the first one is pretty difficult to improve, the second source 
of uncertainty, higher order QED corrections could be well calculated. By these higher 
corrections we mean diagrams presented on Figure 1. Diagrams in the second raw, have 
been calculated in an approximate way in || and are named in the Table I as a muon 
self-energy with electron v. p. Some further theoretical work is necessary for obtaining the 
improved value of the hyperfine splitting of S-levels, since the proton polarizability effect is 
also present there. In summary, having now precise theoretical predictions, the measurement 
of the 2P-2S transition frequency in the muonic hydrogen will lead to the improved value of 
the proton charge radius. 
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TABLES 



correction 


value in meV 


leading order e.v.p. 


205.0074 


rel. corr. to e.v.p . 


0.0594 


double e.v.p. 


0.1509 


two-loop e.v.p. 


1.5079 


three-loop e.v.p. 


0.0076 


muon self-energy + muon v. p. 


-0.6677 


muon self-energy with e.v.p. 


-0.006(1) 


recoil of order a 4 


0.0575 


recoil of order a 5 


-0.0450 


recoil of order a 6 


0.0003 


proton self energy 


-0.0099 


leading finite size of order a 4 


-r 2 5.1974 = -3.862(108) 


finite size of order a 5 


r 3 0.0363 = 0.0232(15) 


finite size of order a 6 


-0.0009(3) 


e.v.p. with finite size 


-r 2 0.0281 = -0.0209(6) 


hadronic v. p. 


0.0113(3) 


proton polarizability 


0.012(2) 


estimate for uncalculated terms 


(0.002) 


sum of corrections to Lamb 


206.085(3) - r 2 5.2255 + r 3 0.0363 


shift in iM with r = 0.862(12) 


= 202.225(108) 



TABLE I. Summary of results for corrections to the Lamb shift 
in muonic hydrogen, e.v.p. denotes the electron vacuum polarization. 



QED corrections are calculated according to Ref. [17 and 
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FIGURES 








FIG. 1. Higher order diagrams contributing to muonic hydrogen Lamb shift. The upper and 
lower horizontal lines denote the muon and the proton respectively. A closed line denotes the 
electron loop. 
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